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ABSTRACT 



The object of this thesis Is to Investigate a number of simple flow 
examples of air-water mixtures and to present the results In a simple quick 
reference for* which can be valuable for engineering use In design and under- 
standing the phenomena of flow of two-component mixtures. The problems 
analysed arei 

1. The one-dimensional trajectory of droplets accelerated In a 
gas stream 

2. The variation of stream properties during the droplet accelera- 
tion process 

3. The variation of stream properties due to well friction with 
water present In the stream. 

The procedure used was entirely analytical. The results of the first 
problem were calculated, tabulated, and plotted in dimensionless form for 
various water and air velocities, droplet diameters, distance moved down the 
duct, and water and air properties. The plots should be valuable for engi- 
neering use in determining the one-dimensional trajectory of spherioal parti- 
cles in a gas strear. An iteration procedure using the plots can be used where 
air velocity, droplet diameter, or stream properties vary during the acceleration. 

The results of the second and third problem were calculated, tabulated, 
and plotted in dimensionless form for various water to eir velocity ratios, 
water to air mass ratee, air Maeh numbers, and stream property ratios. The 
results should be valuable for engineering use in all regions except near 
choking conditions in the stream. Neatr this condition the air velocity increases 
too rapidly for the drepletp to accelerate with the air. Since this was one 
of the assumptions in the analysis the plots are not correct in this region. 

The analysis of the third problem incorporates in it a pseudo-frictional 
term accounting for the effect of momentum exchange of droplets with the duet 
wall. T he exact nature of this effect is not known, but some correlation with 
experimental data is given. 
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I. IHTRCTUCTIOR 

The object of this thesis is to investigate a number of staple flow 
examples of air-water mixtures. The problems to be analysed aret 

(1) the one-dimensional trajectory of droplets 
aeoelerated In a gas stream, 

(2) the variation of stream properties during the 
acceleration process, and 

(3) the variations of stream properties due to 
wall friction with water present in the stream. 

This investigation is undertaken as part of the Aerothermopresso r 
Project at Massachusetts Institute of Technology. The Aerothermopressor 
is a device to Increase the stagnation pressure of a gas stream by reducing 
the stagnation temperature through evaporation of water injected into the 
flow. 

In the work on the Aerothermopressor, there has been no quick refer- 
ence guide to aid in design or in understanding the phenomena of flow of 
air-i«ter mixtures. The purpose of this thesis is to analyse these separate 
effects and to present the results in a simple quick reference form which 
ean be valuable for engineering use in problems concerning flow of air- 
water mixtures. 




II. DROPLET TRAJECTORY 



Procedure 

The example studied is the injection of one droplet at a given initial 
velocity into a moving gas stream. The problem is to determine the subse- 
quent droplet velocities as a function of distance moved downstream and the 
properties of the gas and droplet. 

The assumptions made are : 

(1) the gas velocity is constant. 

(2) the droplet is spherical and of constant 
diameter during the acceleration, and 

(3) the table of drag coefficient vs. Reynolds 
number as given in reference (l) and repeated 
in Table I is correct. 

The details of this analysis are given in the appendix. It was necessary to 
perform a graphical integration of a function of Reynolds number. The results 
of one of these integrations is contained in reference (1). The other inte- 
gration was done using the traposoldal rule. The results of both these inte- 
grations are given in Table II. 

Besults 

Figures I and II are non-dimensional plots of droplet velocity vs. dis- 
tance moved downstream ae a function of gas velocity, droplet diameter, and 
the properties of the gas end droplet. This same information is given in 
Table III. Figure I is for droplet velocity less than gas velocity! Figure 
II is for droplet velocity greater than gae velocity. The plots in Figures I 
and II are normalized to specific initial droplet velocities. For droplet 




Telocity less than gas velocity, the initial droplet velocity is assumed zero. 

For droplet velocity greater than gas velocity, the initial droplet velocity 

is assumed twice the air velocity. For an Initial droplet velocity other 

than the normalized conditions, it is necessary to subtract the distance 

read at the initial condition from that of the final condition. This is 

shown in more detail in the appendix. 

Figures I and II are used in the following manner! From given initial 

droplet velocity, gas velocity, droplet diameter, and gas and droplet proper* 
V V d 

ties, calculate and ”&■" • Enter the plot with these quantities to deter- 
a v a V 

mine the Initial point at the intersection of ordinate and the curve 

Vw r f. 

corresponding to ~tf a • Read on the absicoa the initial value T - r* • 'fov 

v a V d w f v 

follow along the same curve to the intersection of that curve and the 

V - f v a 

desired ^ or — - distance from the initial point. Read off the value of 

T f « ® W W y 

T“ r“ which gives the distance from the initial point with the required — K or 
« r w V \ 

the value of ^ at that distance respectively. 

& 

Figures III and IV are plots of distances required to reach specific 



values of ^ • .5».l as a function of droplet diameter, air Mach number, and 
given air properties assuming initial droplet velocity is zero. Figure III 
is for air stagnation temper* ture of 70°F and stagnation pressure of 1 atmos- 
phere. Figure IV is for air stagnation temperature of 15C0°R and stagna- 
tion pressure of 1 atmosphere. 

If the gas velocity, droplet diameter, or the stream and droplet proper- 
ties are not constant, the solution can be approximated by assuming the varia- 
tion to occur in a stepwise fashion and apply the plots successively over 
each step. This can be done in the ease of evaporation in which a curve of 
droplet diameter vs. droplet velocity or distance may be assumed, or in the 




ease of a converging or diverging section in which, because of known area 
changes, a curve of gas velocity vs. distance can be assumed. 

CPIttlMlm 

The results obtained are given In dimensionless form In as simple a 
manner as possible. The plots should be valuable for engineering use In 
determining one-dimensional trajectory of spherical particles In a gas stream. 
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Procedure 

The example studied is the flow of an air-water mixture past two sections 
1 and 2 In the flow. The problem is to determine the stream properties at 
section 2 as a function of the water velocity at 2, the water rate, and the 
stream properties at section 1. 

The assumptions made are: 

(1) Constant area flow 

(2) Adlabatie— no change in stagnation temperature 

(3) No evaporation of water droplets 

(4) No change in t«mpereture of droplets 

(5) Perfect gae relation holds for air. 

The details of this analysis are given in the appendix. The general method 
of performing the analysis is given in reference (2). 



Results 

T t ® P V 

Values of , py» and —y were calculated for y ■ 0, 1 for 

00 & 

various wster rates and Mach numbers* The star condition is a normalised 

V 

condition corresponding to y ■ 0 and M - 1 end does not depend on the water 

a 

rate. These values are tabulated in Tables IV through VIII. The guiding prin- 

olple in using the tablee le that ijr? is constant in going between tables. 

o 

Figures V through X show the effect of droplet acceleration on stream 

v w V 

properties assuming that at section 1, ~ * 0, and at section 2, ■ 1 for 

a a 

various initial Mach numbers and water rates. The numbers on each curve indi- 

cate the water rate to air rate ratio, — 11 , and the letter S on some branches 

a 

of the curves indicate that a noraal shock in the flow Is necessary to reach 
these states. 
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Figurea V through X are used in the fo' loving manner* Knowing the 

stream properties at section 1, the point of injection of the water, enter 

v 

the required plot with the ^ and ^ . Read on the ordinate the value of the 

A 

ratio of the desired property at section 2 to that at section 1, 

Conclusion 

The plots given are easy to use and are correct within the assumptions. 
However, near choking conditions the air velocity increases very rapidly. It 
is erroneous to assume that the droplets will be able to accelerate with the 
stream in the absenoe of very large drag coefficients. Therefore, near 
choking conditions the plots will not give a oorrect answer to the variation 
of stream properties due to droplet acceleration. However, the value of 
where choking of the flow at section 2 is indicated is approximately correct. 
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tL ^irtEP FAWHD link Amraia 



The example studied Is the flow of an air-water mixture in & pipe with 
wall friction and droplets striking the well. The problem Is, knowing the 
initial stream properties and friction factor, to determine the stream proper- 
ties at any section downstream. 

The assumptions made are i 

(1) Constant area flow 

(2) Adlabatio— no change in stagnation temperature 

(3) Vo evaporation of water droplets 

(A) Wo change in temperature of droplets 

(5) Perfect gas relation holds for air 

... V 

(o) For droplet entrained in air stream j ■ 1 
. . V V . 

(7) For droplet on the wall j ^ ■ 0 

(P) In each differential section an amount of water leaves 
the stream and hits the wall. and an equal amount of 
water leaves the wall and is picked up by the stream 
(7) That this amount of water is proportional to the 
length of section 

The details of this analysis are given in the appendix. The general method 
of performing the analysis Is given in reference (2). 

Results 

TUU.. .f L_ • & • £ . k-. b . h — 

calculated for various water rates and Mach numbers. The star eondltion is 
a normalised condition corresponding to choking of the flow. Theee values, 



- 1 *- 



to elide rule accuracy, are tabulated in Tables XII through TV and plotted in 
figures XI through XVII. These figures show the effect of friction and 
droplet momentum exchange with the walla on the stream properties. The numbers 
on each curve indicate the water rate to air rate ratio, — * • 



The symbol ^ in the term 



[ 4 H 



kjmx * ncor P orate9 as a pseudo-fric- 



tional tern the effect of mcen^ntum exchange of droplets with the wall. This 

w w 

teisn Is a f>inction of the water rate to air rete ratio, ~ , the duct diameter, 

A 

and the droplet diameter and perhaps other stream properties. The nature 
of this function is not known. 'owever, using the data of reference (3) 
for sand together with the low Mach number analysis described in the appen- 

dix, it appears that K * .4 ^ for 1" diameter pipe and sand diameter of 

A 

200 p and 450 p. It appears that the constant of proportionality is slightly 
smaller for Frailer sand diameters. The effect of duct diameter is no known. 

Figures XI through xVII are used in the following manner i Knowing the 
initial stream properties enter Figure XI with the initial ?ach number, M^, 
and ~ and read off the value OM Then knowing j^4 ^ + kJ and 

the length of duct from the initial section 1 to any other aeotion 2, calculate 

[ 4 p * K j *1,2 * Thon c * loulat « [ 4 § * Jp W 2 “ [ 4 V * K ] L max 1 * [ 4 V * K ] h,2 w 
Enter Figure XI with [^4 § ♦ K~J L MT and ^ and read off Mg. Then entering 

the desired plot of properties with K. , Mg, and , read off the value of 

a 

the normalised property at sections 1 and 2. From this form the ratio of the 
property at 2 to that at 1. Knowing the value of the property at 1, the vlue 
at 2 can be calculated . 

Figure XVII is used to determine the Mach number after a normal shook 

v y V 

knowing the Mach number and just before the shock and assuming ^ • 1 on 

w a « w a 

both sides of the shoch. Fnter the plot with initial K and — ■ to find initial 



point. Follow across horisontally with constant until ranching the sans 
v 

curve. Read the value of M after the shock* 
a 

GgnglMAgD 

Provided more information can be obtained concerning the pesudo-frio- 
tlonal tern* K, the plots given can be valuable for engineering use* This 
could well form the basis of a thesis to determine whether this simple model 
of the flow is correct and to correlate data in the field to two-component 
flows. 

In addition* very near choking conditions this model of the flow is not 
oorreot since the air velocity increases very rapidly in a very short length. 
It is erroneous to assume that the droplets can accelerate with the stream 
in the absence of very large drag coefficients. Therefore, near ohoking 
conditions the plots will not give a correct answer to the variation of 
stream properties due to friction* 
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A = CROSS SECTIONAL AREA 
C = SPEEX5 OF SOUND 
Cd = DRAG COEFFICIENT 

Cp= specific neat at CONSTANT pressure: 

D* HYDRAULIC DIAMETER OF DUCT 
dur= DROPLET DIAMETER 
1= FRICTION COEFFICIENT OF DUCT 
F * IMPULSE FUNCTION (PA+FAY 2 ) 
h - SPECIFIC ENTHALPY 
K = RATIO OF SPECIFIC HEATS (<tyCv) 

Lmm- maximum length of duct for CHOKING flow 
tA • MACH NUMBER 
m = MASS OF DROPLET 
P = STATIC PRESSURE 
ft'P»* ISENTROPIC STAGNATION PRESSURE 
R = REYNOLDS NUMBER 
T- TEMPERATURE ( ABSOLUTE ) 

ISENTROPIC STAGNATION TEMPERATURE (ABSOLUTE) 
t = TIME 
V - VELOCITY 
ur= MASS RATE OF FLOW 
X * DISTANCE ALONG DUCT 
f= MASS DENSITY 
l> * KINEMATIC VISCOSITY 

( )*= Refers to a normalized choking condition 

( REFERS TO AIR 
( )ur= REFERS TO WATER 

( )o * REFERS TO ISENTROPIC STAGNATION CONDITION 
( )i = REFERS TO SECTION I 
()i= REFERS TO SECTION Z 
( )(?» REFERS TO CONDITION RELATIVE TO AIR 




DROPLET OF ANALYSIS 

INOlVIOUfiL DR0PLCT5 INJECTED INTO M0V1N6 AIR STREAM 



Let : \/a » Constant 

DeFiNC: l/fl * \J a - \/ur (TOR l/a, ^ l/af ) 

\ 4 r - \A*- (> 0 /? lk< l/urjl 

D _ Aft <^ur \Z<R 
~ //« 

ro/?cc eooMio^ 0 n DROPier 

M ^ s - Co A \/p 
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lEcfur*) 



dt 



h* %TT<iv 
Ml * g" n dur 





DROPLET TRAJECTORY - CONT 




INTEAMS OF REAL DISTANCES 

X .,2 = V a t h2 +- X „,. 4 (3) 

NOTE: IF 14 > Yu - , X* A5 Af Mt&AT/l/E DtREC TICN 
IN FO Q AT ICV (4) SI&N OF LAST TEW IS: 



H- IF Va < Vur 

1 T 1/a ^ 




Normalizing equation (4) so that initially, 
at point i, i/« = v«. j so that: 

Vur - 0 rO/T 1/a > 14*- 

\/ a , =■ ZV* FOR Vo. < l /««- 

LET THU NORMALIZED POSITION BE POINT O, AND 
ANY OTHER POINT DOWNSTREAM BE POINT X. 




DROPLET TRAJECTORY - CON'T 
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FOR ANY INJECTED V* OTHER THAN THE 
NORMALIZED CONDITION: 

/k _ X c , / 2, fu. X c> y i To. 

i{ us fjr J.ur fv dv T or 



Xo, 



.ur 



A- H/HS e£* 7 V CALCULATED ( TABLE ITT) 

fur 



AND PLOTTED (FIGURES I AMD IT) FOR VALUES OF 

V*. dur from STD ZOOO AND 1*L FROM .01 TO / . 
'io, TZ 



ifll'JLfi X 
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DiscmiiNurn malhis 



EFFECT OF DROPLET ACCELERATION ON AIR STREAM 



ASSUMPTION S ! 

/. ANABATIC - NO CHANGE IN STAGNATION TEMPERATURE 
2. NO EVAPORATION OF WATER DROPLETS 
=3 NO CHANGE IN TEMPERATURE OF WATER DROPLETS 
4. PERFECT GAS RELATION HOLDS FOR AIR 
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MOMENTUM < 



P, * ^ Vc * -S&t Vur, = a * Sf* Viz + ^ Vn 

n * [<'SE fe] «.•»*[<** fey ft. «. 
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NORMALIZING ALL EQUATIONS TO * CONDITION 
WHERE; -ky, - c ; M- I 





ISENTROPIC STAGNATION CONDITIONS' 

I. AS5UM/N6 DROPLETS OECeieRATE- ISEAITROT ICAUY. 
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DISCONTINUITY ANALY 
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2 . ASSUMING DROPLETS ARE NOT DECELERATED 
•Pe = rf'tVM*!* 

. A _ [,+ *£ M 'J 96 * Z X" 

"W “ n Tkf^^JivR [K+i]*-< 

VALUES OF ■&, -& ; , £ , X /»ff£ TABULATED IN 
TABLES JT THR0U6HUJK FOR M FROM O TO Z J 
= . f, . 2, . 5, . ^ | ^ = O /l/W> / ; /i --/ .4 . 

'VALUES OF M , , Jki- , jfk* , _£L Tz I/a 

* ’ p°, ‘ > p, > > v, 

ARE PLOTTED VS Mi I IV FIGURES X THROUGH X 
4/VD /W 2 /5 TABULATED KS /W, //»/ TABLE LF FOR 
M, FROM OTO Z j //=/.+; 

ASSUMING NgL = O /WD = / . 
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MODIFIED FAN NO LINE ANALYSIS 

FRICTIONAL, ADIABATIC , CONSTANT-AREA . 

ASSUMPTIONS : 

I. ADIABATIC- No CHANGE IN STAGNATION TEMPERATURE 
Z.NO evaporation OF water DROPLETS 

3. NO CHANGE IN TEMPERATURE OF WATER DROPLETS 

4. PERFECT GAG RELATION HOLDS FOR AIR 

5. FOR DROPLET ENTRAINED IN AIR STREAM; 

G. FOR DROPLET ON THE WALL ; - © 

TIN EACH DIFFERENTIAL SECTION, LENGTH dt , An 

amount of water, dw, , leaves the stream 
and purs the wall, in addition an amount of 

water , dur z , LEAVES THE WALL AND IS PICKED 
UP BY THE STREAM. 



CONTINUITY-. 



/>!/* CONSTANT 

dl + 

f V 

EQUATION OF STATE: 

P= f RT 

dP _ d_£- + dT~ 
P f ~T 



definition of m ■■ . 

£ - 



Kft T 

dM z _ dV z _ dr 

M z \l z T 





(Jj 



MODIFIED FAN NO LINE ANALYSIS - COM'r 



ENERGY-. 




uK.Cha.i-dht) 

(UTatUTur -d&i + dwz) 

( y % + d VX) 



O = Ufa. dha. + (u/a. 1 -Ufw) dlC’ - ( dor, -dur i )(M?' V- AM 

DROPPING ZND ORDER DIFFERENTIAL TERM AND 
REARRANGING) NOTING THAT: 



dha * Cpo. dT j 

Ml = m * m 

Z 2. > 

HR - Ji± 

2C? 2 < 

- dur^ = duTur 



0 




&(»&)** dot. 

t+¥U*W& M * + 



/V STEADY state \NITH WATER on THE WALLS, 
ASSUME THE FILM THICKNESS IS CONSTANT WITH 
TIME AT ANY SECTION AND THE FILM \IELOCITY 
IS NE6LI6IBLE . 



MODIFIED FANNO LINE ANALYSIS- CON'T 



MOMENTUMj 
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To rTTDdY dur,)l d»> *(0) 



-AdP-TrTTDe/K = (wi.TurJ)dV + dur t V~(dur, -dur^dV 




DEFINE: 

r Tor _ _TV_ 

T " -fc /V* fc KPM* 




ASStfMf TtfAf //V AW SPECIFIC CASE THE AMOUNT OF 
W ATE A PICKED UP BY THE STREAM FROM THE I NALLS 
IS PROPORTIONAL TO THE LENGTH. 

2 dJE*. * K dX 



RESULTS OF SIMILTAHEOUS SOLUTION OF EQUATIONS 
(I) THROUGH (S) ARE PRESENTED IN TABLE X 



MODIFIED FANNO LINE ANAL Y5IS- CON'T 

TABLE X 
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results of integration and integral relations 

ARE GIVEN IN TABLE IXT 










MODIFIED FANNO LINE ANALYSIS - CON'T 

TABLE JL 



FORMULAS IN INTEGRAL FORM 




ASSUMING DROPLETS DO NOT DECELERATE: 





mpiEM fan no line analysis - con't 
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